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Millimeter Wave MIMO with Lens Antenna Array: 
A New Path Division Multiplexing Paradigm 

Yong Zeng and Rui Zhang 


Abstract —Millimeter wave (mmWave) communication over the 
largely unused mmWave spectrum is a promising technology 
for the fifth-generation (5G) cellular systems. To compensate 
for the severe path loss in mmWave communications, large 
antenna arrays are generally used at both the transmitter and 
receiver to achieve significant heamforming gains. However, 
the high hardware and power consumption cost due to the 
large number of radio frequency (RF) chains required renders 
the traditional beamforming method impractical for mmWave 
systems. It is thus practically valuable to achieve the large- 
antenna gains, but with only limited number of RF chains 
for mmWave communications. To this end, we study in this 
paper a new lens antenna array enabled mmWave multiple-input 
multiple-output (MIMO) communication system. We first show 
that the array response of the proposed lens antenna array at the 
receiver/transmitter follows a “sine” function, where the antenna 
with the peak response is determined by the angle of arrival 
(AoA)/departure (AoD) of the received/transmitted signal. By 
exploiting this unique property of lens antenna arrays along with 
the multi-path sparsity of mmWave channels, we propose a novel 
low-cost and capacity-achieving MIMO transmission scheme, 
termed orthogonal path division multiplexing (OPDM). With 
OPDM, multiple data streams are simultaneously transmitted 
in parallel over different propagation paths with simple per- 
path signal processing at both the transmitter and receiver. 
For channels with insufficiently separated AoAs and/or AoDs, 
we also propose a simple path grouping technique with group- 
based small-scale MIMO processing to mitigate the inter-path 
interference. Numerical results are provided to compare the 
performance of the proposed lens antenna arrays for mmWave 
MIMO system against that of conventional arrays, under differ¬ 
ent practical setups. It is shown that the proposed system achieves 
significant throughput gain as well as complexity and hardware 
cost reduction, both making it an appealing new paradigm for 
mmWave MIMO communications. 

Index Terms —Lens antenna array, millimeter wave communi¬ 
cation, antenna selection, path division multiplexing, inter-path 
interference. 


I. Introduction 

The fifth-generation (5G) wireless communication systems 
on the roadmap are expected to provide at least 1000 times ca¬ 
pacity increase over the current 4G systems ||T|. To achieve this 
goal, various technologies have been proposed and extensively 
investigated during the past few years Among others, 
wireless communication over the largely unused millimeter 
wave (mmWave) spectrum (say, 30-300GHz) is regarded as 
a key enabling technology for 5G and has drawn significant 
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interests recently (see Q-fhl and the references therein). Ex¬ 
isting mmWave communication systems are designed mainly 
for short-range line-of-sight (LOS) indoor applications, e.g., 
wireless personal area networking (WPAN) Q and wireless 
local area networking (WLAN) Q. While recent measurement 
results have shown that, even in non-line-of-sight (NLOS) out¬ 
door environment, mmWave signals with satisfactory strengths 
can be received up to 200 meters Q, pO) , which indicates 
that mmWave communications may also be feasible for future 
cellular networks with relatively small cell coverage. 

MmWave signals generally experience orders-of-magnitude 
more path loss than those at much lower frequency in 
existing cellular systems. On the other hand, their smaller 
wavelengths make it practically feasible to pack a large 
number of antennas with reasonable form factors at both the 
transmitter and receiver, whereby efficient MIMO (multiple- 
input multiple-output) beamforming techniques can be applied 
to achieve highly directional communication to compensate for 
the severe path loss pOj- GD- However, traditional MIMO 
beamforming is usually implemented digitally at baseband 
and thus requires one dedicated radio frequency (RL) chain 
for each transmit/receive antenna, which may not be feasible 
in mmWave systems due to the high hardware and power 
consumption cost of the large number of RL chains required. 
To reduce the cost and yet achieve the high array gain, analog 
beamforming has been proposed for mmWave communications 
Gg-ig, which can be implemented via phase shifters in the 
RL frontend, and thus requires only one RL chain for the entire 
transmitter/receiver. Despite of the notable cost reduction, 
analog beamforming usually incurs significant performance 
loss due to the constant-amplitude beamformer constraint 
imposed by the phase shifters, as well as its inability to 
perform spatial multiplexing for high-rate transmission. To 
enable spatial multiplexing, hybrid analog/digital precoding 
has been recently proposed l[T3-@, where the precoding 
is implemented in two stages with a baseband digital pre¬ 
coding using a limited number of RL chains followed by 
a RL-band analog processing through a network of phase 
shifters. Since the hybrid precoding in general requires a large 
number of phase shifters, antenna subset selection has been 
proposed in [ |23| by replacing the phase shifters with switches. 
However, antenna selection may cause significant performance 
degradation due to the limited array gains resulted p4| , p5) , 
especially in highly correlated MIMO channels as in mmWave 
systems. 

Besides, another promising line of research for mmWave 
or large MIMO systems aims to reduce signal processing 
complexity and RL chain cost without notable performance 
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degradation by utilizing advanced antenna designs, such as 
the lens antenna array p^-|[30l. As shown in Fig. [T] a 
lens antenna array is in general composed of two main 
components; an electromagnetic (EM) lens and a matching 
antenna array with elements located in the focal region of 
the lens. Generally speaking, EM lenses can be implemented 
via three main technologies: i) the dielectric lenses made of 
dielectric materials with carefully designed front and/or rear 
surfaces |3T), p2| ; ii) the traditional planar lenses consisting 
of arrays of transmitting and receiving antennas connected via 
transmission lines with variable lengths p^ , |j^; and iii) the 
modern planar lenses composed of sub-wavelength periodic 
inductive and capacitive structures p5) , p^ . Regardless of 
the actual implementation methods, the fundamental principle 
of EM lenses is to provide variable phase shifting for EM rays 
at different points on the lens aperture so as to achieve angle 
of arrival (AoA)/departure (AoD)-dependent energy focusing, 
i.e., a receiving (transmitting) lens antenna array is able to 
focus (steer) the incident (departure) signals with sufficiently 
separated AoAs (AoDs) to (from) different antenna subsets. 
In p7) , the concept of beamspace MIMO communication 
is introduced, where the lens antenna arrays are used to 
approximately transform the signals in antenna space to 
beamspace, which has much lower dimensions, to significantly 
reduce the number of RE chains required. However, the studies 
in HZ) focus on the LOS mmWave channels, where spatial 
multiplexing is possible only for very short transmission range 
(e.g. a few meters) and/or extremely large antenna apertures. 
In a parallel work p^ , the lens antenna array is applied 
to the massive MIMO cellular system with large number 
of antennas at the base station (BS) p7|-p9), which is 
shown to achieve significant performance gains as well as cost 
reduction as compared to the conventional arrays without lens. 
However, the result in p9) is only applicable for the single¬ 
input multiple-output (SIMO) uplink transmission, instead of 
the more general setup with lens antenna arrays applied at both 
the transmitter and receiver. Moreover, neither iz) nor | |29l 
fully explores the characteristics of mmWave channels, such 
as the multi-path sparsity Q due to limited scattering and the 
frequency selectivity in broadband transmission. 

In this paper, we study the mmWave MIMO communication 
where both the transmitter and receiver are equipped with 
lens antenna arrays. Due to the AoA/AoD-dependent energy 
focusing, in mmWave systems with limited number of multi¬ 
paths, the signal power is generally focused on only a small 
subset of the antenna elements in the lens array; as a result, 
antenna selection can be applied to significantly reduce the 
RE chain cost, yet without notably comprising the system 
performance, which is in sharp contrast to the case of applying 
antenna selection with the conventional arrays p4) , p5) . 
Eurthermore, for mmWave channels with sufficiently separated 
AoAs/AoDs, different signal paths can be differentiated in 
the spatial domain with the use of the lens antenna array. 
Therefore, the detrimental multi-path effect in wide-band 
communications, i.e., the inter-symbol interference (ISI), can 
be easily alleviated in the lens array MIMO systems, without 
the need of sophisticated ISI mitigation techniques such as 
equalization, spread spectrum, or multi-carrier transmission 


pBl . In fact, in the favorable scenario where the AoAs/AoDs 
are sufficiently separated, the lens MIMO system can be 
shown to be equivalent to a set of parallel additive white 
Gaussian noise (AWGN) sub-channels, each corresponding 
to one of the multi-paths, for both narrow-band and wide¬ 
band communications. Thus, multiple data streams can be 
simultaneously multiplexed and transmitted over these sub¬ 
channels in parallel, each over one of the multi-paths with 
simple per-path processing. We term this new MIMO spatial 
multiplexing scheme enabled by the lens antenna array as 
orthogonal path division multiplexing (OPDM), in contrast to 
the conventional multiplexing techniques over orthogonal time 
or frequency)^ We summarize the main contributions of this 
paper as follows. 

• Eirst, we present the array configuration for the proposed 
lens antenna array in detail, and derive its corresponding 
array response. Our result shows that, different from the 
conventional arrays whose response is generally given 
by phase shifting across the antenna elements, the array 
response for the lens antennas follows a “sine” function, 
where the antenna with peak response is determined by 
the AoA/AoD of the received/transmitted signal. This 
analytical result is consistent with that reported in prior 
works based on simulations |30|, |[32) or experiments 


j. With the derived array response, the channel model 
for the lens MIMO system is obtained, which is compared 
with that of a benchmark system using the conventional 
uniform planar arrays (UPAs). 

• Next, to obtain fundamental limit and draw insight, we 
consider the so-called “ideal” AoA/AoD environment, 
where the signal power of each multi-path is focused 
on one single element of the lens array at the re¬ 
ceiver/transmitter. We show that the channel capacity in 
this case is achieved by the novel OPDM scheme, which 
can be easily implemented by antenna selection with 
only L transmitting/receiving RE chains, with L denoting 
the number of multi-paths. Notice that L is usually 
much smaller than the number of transmitting/receiving 
antennas in mmWave MIMO channels due to the multi- 
path sparsity. We further compare the lens array based 
mmWave MIMO system with that based on the conven¬ 
tional UPAs, in terms of capacity performance as well as 
signal processing complexity and RE chain cost. 

• Einally, the mmWave lens MIMO is studied under the 
practical setup with multi-paths of arbitrary AoAs/AoDs. 
We propose a low-complexity transceiver design based 
on path-division multiplexing (PDM), applicable for both 
narrow-band and wide-band communications, with per- 
path maximal ratio transmission (MRT) at the transmitter 
and maximal ratio combining (MRC)/minimum mean 
square error (MMSE) beamforming at the receiver. We 

• Note that OPDM also differs from the conventional sectorized antenna and 
space-division-multiple-access (SDMA) techniques. Although they similarly 
exploit the different AoAs/AoDs of multiuser/multi-path signals, the former 
achieves spatial signal separation only in a coarse scale (say, 120 degrees 
with a 3-sector antenna anay), while the latter obtains finer spatial resolution 
but with sophisticated beamforming/precoding. In contrast, with the proposed 
OPDM, high spatial resolution is achieved without the need of complex array 
signal processing. 




3 


analytically show that in the case of wide-hand com¬ 
munications, the proposed design achieves perfect ISl 
rejection if either the AoAs or AoDs (not necessarily 
both) of the multi-path signals are sufficiently separated, 
which usually holds in practice. Moreover, for cases 
with insufficiently separated AoAs and/or AoDs, we 
propose a simple path grouping technique with group- 
based small-scale MIMO processing to mitigate the inter¬ 
path interference. 

It is worth pointing out that there has been an upsurge of 
interest recently in exploiting the angular domain of multi¬ 
path/multiuser signals in the design of massive MIMO sys¬ 
tems. For example, by utilizing the fact that there is limited an¬ 
gular spread for signals sent from the mobile users, the authors 
in EH propose a channel covariance-based pilot assignment 
strategy to mitigate the pilot contamination problem in multi¬ 
cell massive MIMO systems. Similarly in [ |42) , | |43] , an AoA- 
based user grouping technique is proposed, which leads to the 
so-called joint spatial division and multiplexing scheme that 
makes massive MIMO also possible for frequency division 
duplexing (FDD) systems due to the significantly reduced 
channel estimation overhead after user grouping. In | [44j , an 
OFDM (orthogonal frequency division multiplexing) based 
beam division multiple access scheme is proposed for massive 
MIMO systems by simultaneously serving users with different 
beams at each frequency sub-channel. In this paper, we also 
exploit the different AoAs/AoDs of multi-path signals for 
complexity and cost reduction in mmWave MIMO systems, by 
utilizing the novel lens antenna arrays at both the transmitter 
and receiver. 

The rest of this paper is organized as follows. Section |II] 
presents the array architecture as well as the array response 
function of the proposed lens antenna, based on which the 
MIMO channel model for mmWave communications is de¬ 
rived. The benchmark system using the conventional UPAs is 
also presented. In Section m we consider the case of “ideal” 
AoA/AoD environment to introduce OPDM and demonstrate 
the great advantages of applying lens antenna arrays over con¬ 
ventional UPAs in mmWave communications. In Section |IV] 
the practical scenario with arbitrary AoAs/AoDs is considered, 
where a simple transceiver design termed PDM applicable 
for both narrow-band and wide-hand communications is pre¬ 
sented, and a path grouping technique is proposed to further 
improve the performance. Finally, we conclude the paper and 
point out future research directions in Section |V] 

Notations: In this paper, scalars are denoted by italic letters. 
Boldface lower- and upper-case letters denote vectors and 
matrices, respectively. denotes the space of M x N 

complex-valued matrices, and I represents an identity matrix. 
For an arbitrary-size matrix A, its complex conjugate, trans¬ 
pose, and Hermitian transpose are denoted by A*, A^, and 
A^, respectively. For a vector a, ||a|| denotes its Euclidean 
norm, and diag(a) represents a diagonal matrix with the 
diagonal elements given in a. For a non-singular square 
matrix S, its matrix inverse is denoted as S“^. The symbol 
j represents the imaginary unit of complex numbers, with 
= —1. The notation ® denotes the linear convolution 
operation. i5(-) denotes the Dirac delta function, and sinc(-) 


is the “sine” function defined as sinc(a;) = sin(7rx)/(7ra;). For 
a real number a, [aj denotes the largest integer no greater 
than a, and round (a) represents the nearest integer of a. 
Furthermore, U[a,b] represents the uniform distribution in 
the interval [a,b]. Af{p,C) and CAf{fi,C) denote the real¬ 
valued Gaussian and the circularly symmetric complex-valued 
Gaussian (CSCG) distributions with mean p and covariance 
matrix C, respectively. For a set S, |5| denotes its cardinality. 
Furthermore, n S 2 and U S 2 denote the intersection and 
union of sets Si and S 2 , respectively. 

II. System Description and Channel Model 
A. Lens Antenna Array 

A lens antenna array in general consists of an EM lens and 
an antenna array with elements located in the focal region of 
the lens. Without loss of generality, we assume that a planar 
EM lens with negligible thickness and of size Dy x is 
placed on the y-z plane and centered at the origin, as shown 
in Eig. By considering only the azimuth AoAs and AoDs0 
the array elements are assumed to be placed on the focal 
arc of the lens, which is defined as a semi-circle around the 
lens’s center in the azimuth plane (i.e., x-y plane shown in 
Fig-B with radius F, where F is known as the focal length of 
the lens. Therefore, the antenna locations relative to the lens 
center can be parameterized as Bm{xm = F cos 9m, Um = 
—F srnOmyZm = 0), where 9m G [—7r/2,7r/2] is the angle of 
the mth antenna element relative to the x-axis, m G A4, with 
M = {0,±1, • • • ,±(M — l)/2} denoting the set of antenna 
indices and M representing the total number of antennas. 
Note that we have assumed that M is an odd number for 
convenience. Eurthermore, we assume the so-called critical 
antenna spacing, i.e., the antenna elements are deployed on 
the focal arc so that {9m — sin^m} are equally spaced in the 
interval [—1,1] as 

" TTl 

9m = mG M, (1) 

where D = Dy/X is the effective lens dimension along the az¬ 
imuth plane, with A denoting the carrier wavelength. It follows 
from Q that M and D are related via M = 1 -f [2Z)j, i.e., 
more antennas should be deployed for larger lens dimension 
D. It is worth mentioning that with the array configuration 
specihed in Q, antennas are more densely deployed in the 
center of the array than those on each of the two edges. 

We first study the receive array response by assuming that 
the lens antenna array is illuminated by a uniform plane 
wave with AoA as shown in Eig. B Denote by xo{4>) the 
impinging signal at the reference point (say, the lens center) 
on the lens aperture, and Cmif) the resulting signal received 
by the mth element of the antenna array, m G Ai. The array 
response vector a.{(j)) G whose elements are defined 

by the ratio am(f) = 't’m{4>)/xo{(j)), can then be obtained in 
the following lemma. 

^For simplicity, we assume that the elevation AoAs/AoDs are all zeros, 
which is practically valid if the height difference between the transmitter and 
the receiver is much smaller than their separation distance. 
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Fig. 1: The schematic diagram of a lens antenna array with 
an incident uniform plane wave with AoA (j). 


Lemma 1: For the lens antenna array with critical antenna 
spacing as specified in Q, the receive array response vector 
a.{(j)) as a function of the AoA (j) is given by 

am{(j>) = '/Asmc{m — D^), m G M, (2) 

where A = DyD^/X^ is the effective aperture of the EM lens, 
and (j) = sini^ G [—1,1] is referred to as the spatial frequency 
corresponding to the AoA f. 

Proof: Please refer to Appendix ■ 

Different from the traditional antenna arrays without lens, 
whose array responses are generally given by the simple phase 
shifting across different antenna elements (see e.g. ( [TT] l for 
the case of UPAs), the “sinc”-function array response in (|^ 
demonstrates the AoA-dependent energy-focusing capability 
of the lens antenna arrays, which is illustrated in Fig. 
Specifically, for any incident signal with a given AoA f, the 
received power is magnified by approximately A times for 
the receiving antenna located in the close vicinity of the focal 
point D(j)\ whereas it is almost negligible for those antennas 
located far away from the focal point, i.e., antennas with 
\m — D(j>\ ^ 1. As a result, any two simultaneously received 
signals with sufficiently different AoAs f and f such that 
\4> ~ 4''\ ^ 1/-0 can be effectively separated in the spatial 
domain, as illustrated in Fig. [^assuming a lens antenna array 
with A = 100 and D — IQ for two AoAs with sinf — 0 
and 0.18, respectively. Thus, we term the quantity l/D as the 
array’s spatial frequency resolution, or approximately the AoA 
resolution for large D 

On the other hand, since the EM lens is a passive device, 
reciprocity holds between the incoming and outgoing signals 
through it. As a result, the transmit response vector for steering 
a signal towards the AoD f can be similarly obtained by 
Lemma [T] 

B. Channel Model for MmWave Lens MIMO 

In this subsection, we present the channel model for the 
mmWave lens MIMO system, where both the transmitter and 



Eig. 2: Array response of a lens antenna array with A = 100 
and D = IQ for two different AoAs. 



Eig. 3: A mmWave lens MIMO system in multi-path environ¬ 
ment. 


receiver are equipped with lens antenna arrays with Q and M 
elements, respectively, as shown in Eig. Under the general 
multi-path environment, the channel impulse response can be 
modeled as 

L 

H(f) = ^ aia.R{(j)R^i)a^- n), (3) 

;=i 


where H(f) is an MxQ matrix with elements hmq{t) denoting 
the channel impulse response from transmitting antenna q G Q 
to receiving antenna m G A4, with Q and A4 respectively 
denoting the sets of the transmitting and receiving antenna 
indices as similarly defined in Section |II-A L denotes the 
total number of significant multi-paths, which is usually small 
due to the multi-path sparsity in mmWave channels 0 ; cei 
and Ti denote the complex-valued path gain and the delay 
for the Zth path, respectively; and are the azimuth 
AoA and AoD for path I, respectively; and G and 

ax G represent the array response vectors for the lens 

antenna arrays at the receiver and the transmitter, respectively. 
Note that in ([^, we have assumed that the distances between 
the scatterers and the transmitter/receiver are much larger than 
the array dimensions, so that each multi-path signal can be 
well approximated as a uniform plane wave. 


' i / 
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Denote by At and An the effective lens apertures, and Dt 
and Dn the lens’s effective azimuth dimensions at the trans¬ 
mitter and at the receiver, respectively. Based on Lemma 
the elements in the receive and transmit array response vectors 
Sin and can be respectively expressed as 


aR,m{ 4 >R,i) = \A 4 ^sinc(m - mG M, ( 4 ) 

aT,q{ 4 >T,i) = \/^sinc(g — Dt^t,i), q G Q, (5) 

where (j)R^i = sm{cl)R^i) and = sm{(j)T,i) are the 

AoA/AoD spatial frequencies of the Zth path. Without loss 
of generality, (j)T,i G [—1,1] of the L multi-paths can 
be expressed in terms of the spatial frequency resolutions 
associated with the receiving/transmitting arrays as 


~ mi + crj 2 

9R,i = -tx-) 9T,l 

Ur 


qi + eT,i 
Dt 


,L, 


( 6 ) 


where mi G A4 and qi G Q are integers given by m; = 
round((^ij^;i)i^) and qi = roundand cr^ and ct,; 
are fractional numbers in the interval [—1/2,1/2]. Intuitively, 
mi (or qi) in gives the receiving (transmitting) antenna 
index that is nearest to the focusing point corresponding to the 
AoA (AoD) of the hh path; whereas ; and et,; represent 
the misalignment from the exact focusing point of the Zth path 
signal relative to its nearest receiving/transmitting antenna. By 
substituting (|^ into Q and (|^, the channel impulse response 
in © can be equivalently expressed as 


L 

hmq{t) ='^ai\/A~i^^smc{m-mi - erj) 

1=1 ’ 
X smc{q - qi- €T,i)S{t - n), m G M,q G Q. 

Loosely speaking, 0 implies that the signal sent by the 
transmitting antenna with index q = qi will be directed towards 
the receiver mainly along the Zth path, and be mainly focused 
on the receiving antenna with index m = mi, as illustrated in 
Fig.[3l 

With the channel impulse response matrix H(f) given in 
0, the baseband equivalent signal received by the receiving 
lens antenna array can be expressed as 


r{t) = H(f) ® x(f) -b z{t) 

L 

= - n) + z(i), (8) 

where x(f) G denotes the signal sent from the Q 

transmitting antennas, and z{t) G represents the 

AWGN vector at the receiving antenna array. In the special 
case of narrow-band communications where the maximum 
excessive delay of the multi-path signals is much smaller than 
the symbol duration Tg, i.e., maxjri — r//| <C Tg « XjW 

with W denoting the signal bandwidth, we have ti k, t and 
x{t — Ti) « x(< —r), yi. As a result, by assuming perfect time 
synchronization at the receiver, the general signal model for 
the wide-band communications in 0 reduces to 

r(f) = Hx(f)-bz(f), (9) 

where H = denotes the narrow- 

band MIMO channel. 




i]1 '4 /f ’ / 
4 / t / V' 

Y I ^ U 
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(b) Uniform planar array 

Fig. 4: 3D schematic diagrams of a lens antenna array versus 
an UFA with the same physical dimensions. 


C. Benchmark System: MmWave MIMO with Uniform Planar 
Array 

As a benchmark system for comparison, we consider the 
mmWave communications in the traditional MIMO setup 
employing conventional antenna arrays without the EM lens. 
In particular, we assume that the transmitter and the receiver 
are both equipped with the UPAs with Qu and Mu elements, 
respectively, with adjacent elements separated by distance 
di) = 0.5A. For fair comparison, we assume that Qu and 
Mu are designed such that the UFA has the same physical 
dimensions (or equivalently the same effective apertures At 
and Ar) as the lens antenna of interest, as illustrated in Fig. 
Accordingly, it can be obtained that Qu = DyDz/du = 
4 At > Q and Mu — 4 Ar > M, i.e., in general more 
antennas need to be deployed in the conventional UFA than 
that in the lens antenna array to achieve the same array 
aperture, since the energy focusing capability of the EM lens 
effectively reduces the number of antenna elements required 
in lens array. This may compensate the additional cost of 
EM lens production and integration in practice. Denote by 
Hu(f) G the channel impulse response matrix in the 












6 


mmWave MIMO with UPAs. We then have 

L 

Hu(i) = X! - Ti), (10) 

1=1 

where and (j)T,i are dehned in (|^, and and 

are the array response vectors corresponding to the UPAs 
at the receiver and transmitter, respectively, which are given 
by phase shifting across different antenna elements as p3] 


^R,vi4’) = 

1 

V 

1 , 


gi^ML 

i<t>) 

T 

5 

( 11 ) 

aT,u((/') = 

/ -^T 

V Qu 

1, 

gi^2(0)^ ... 

gi't’Qu 

w' 

T 

•) 

( 12 ) 

with 1X1 

= 2 ,.. 

,Mv or 2 ,- 

■ ,Qv, 

denoting 

the 


phase shift of the mth array element relative to the hrst 
antenna. The input-output relationships for the UPA-based 
wide-band/narrow-band mmWave MIMO communications can 
be similarly obtained as in (|^ and 0 , respectively, and are 
thus omitted for brevity. 

In this paper, we assume that the MIMO channel is perfectly 
known at the transmitter and receiver for both the proposed 
lens MIMO and the benchmark UPA-based MIMO systems. 

III. Lens MIMO under Ideal AoAs and AoDs 

To demonstrate the fundamental gains of the lens MIMO 
based mmWave communication, we first consider an “ideal” 
multi-path propagation environment, where the spatial fre¬ 
quencies corresponding to the AoAs/AoDs 

of the L paths are all integer multiples of the spatial fre¬ 
quency resolutions of the receiving/transmitting lenses, i.e., 
{^R,i,^T,i}iLi dehned in ([^ are all zeros. Furthermore, we 
assume that all the L signal paths have distinct AoAs/AoDs 
such that mi' ^ m./ and qi' ^ qi, VT 7 ^ 1. In this case, 
we show that the multi-path signals in the lens antenna 
enabled mmWave MIMO system can be perfectly resolved 
in the spatial domain, thus leading to a new and capacity- 
achieving spatial multiplexing technique called OPDM. We 
also show that with OPDM, the lens antenna based mmWave 
MIMO system achieves the same (or even better) capacity 
performance in the narrow-band (wide-hand) communications 
as compared to the conventional UPA based mmWave MIMO, 
but with dramatically reduced signal processing complexity 
and RF chain cost. 

A. Orthogonal Path Division Multiplexing 

In the “ideal” AoA/AoD environment as dehned above, the 
channel impulse response from the transmitting antenna q to 
receiving antenna m given in o reduces to 

L 

hmq{t) = X! oqyj ARAT5{m - mi)6{q - qi)S{t - n). 

1=1 

(13) 

The expression in implies that the signal transmitted 
by antenna q will be received at antenna m if and only if 
there exists a propagation path such that the focusing points 
corresponding to its AoA and AoD align exactly with the 
locations of antenna m and q, respectively, i.e., m = mi and 




Fig. 5: Equivalent input-output relationship for OPDM. 

q = qi. Denote by Xq{t) the signal sent by antenna q of the 
transmitting lens array, where q G Q. The signal received by 
antenna m (by ignoring additive noise for the time being) can 
then be expressed as 

r-mit) = ^ hmq{t) @ Xq{t) 

qeQ 

AfiAxaiXq^ {t — Ti), if m = mi for some I, 

0 , otherwise. 

(14) 

Under the assumption of perfect time synchronization at each 
of the receiving antennas, i.e., r/ is known at the receiver and 
perfectly compensated at antenna mi, ( [T4l i can be equivalently 
written as 

Tmi \/Aj^Aj'dlX -f , I 1, * * • , L, (15) 

where Zmi denotes the AWGN at receiving antenna mi. 
Therefore, the original multi-path MIMO channel has been 
decoupled into L parallel SISO AWGN channels, each cor¬ 
responding to one of the L multi-paths. It is worth men¬ 
tioning that the channel decomposition in ([BJ holds for 
both the narrow-band and wide-hand communications. This 
thus enables a new low-complexity and cost-effective way to 
implement MIMO spatial multiplexing, by multiplexing L data 
streams each over one of the L multi-paths independently, 
which we term as OPDM. 

A schematic diagram of the equivalent input-output rela¬ 
tionship for OPDM is shown in Fig. It is straightforward to 
show that by applying the standard water-hlling (WF) power 
allocation | [40| over each of the L parallel sub-channels with 
power gains {\ai\'^AfiAryjGi, the capacity of the mmWave 
lens MIMO system can be achieved for both narrow-band and 
wide-hand communications. 

B. Capacity Comparison 

Next, we provide capacity comparison by simulations for 
the proposed lens MIMO versus the conventional UPA-based 
MIMO in mmWave communications. For the lens MIMO 
system, we assume that the transmitter and receiver lens 
apertures are both given by Ax = Afi = 20 , and the 
effective azimuth lens dimensions are Dt = Dr = 10 , 
which corresponds to the number of transmitting/receiving 
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antennas as M = Q = 21. For fair comparison, the UPA-based 
MIMO system is assumed to have the same array apertures 
as the lens MIMO, which thus needs Mu = Qu = 80 


transmitting/receiving antennas, as discussed in Section II-C 
We consider a mmWave channel of L = 3 paths, which is 
typical in mmWave communications 0 - We assume a set of 
ideal AoAs/AoDs with (/>tj = G {0, ±0.2}. Furthermore, 
the complex-valued path gains are modeled as ai = 




1 = 1 ,- 




0 , where /3 denotes the large- 


scale attenuation including distance-dependent path loss and 
shadowing, ki represents the power fractional ratio for the 
kh path, with — 1 ’ Vi ^ U[0, 27r] denotes the 

phase shift of the kh path. The value of /3 is set based on 
the generic model —/3dB = Ci ± 10 c 2 logiQ(d) + where 
Cl and C 2 are the model parameters, d is the communication 
distance in meters, and ^ ~ A/^( 0 , e^) denotes the lognormal 
shadowing. We assume that the system is operated at the 
mmWave frequency / = 73 GHz, for which extensive channel 
measurements have been performed and the model parameters 
have been obtained as ci = 86 . 6 , C 2 = 2.45, and e = 8 dB 0. 
Furthermore, we assume d — 100 meters, with which the path 
loss is 136 dB, or E[/3] = —136 dB, with the expectation taken 
over the log-normal shadowing. In addition, the multi-path 

can be modeled as ki = , 

-ilQ-o.iz^^ where Uk ~ t/[ 0 , 1 ] and z\ ~ 
A/^( 0 ,(^^) are random variables accounting for the variations 
in delay and in lognormal shadowing among different paths, 
respectively 0. For mmWave channels at / = 73 GHz, r^. and 
C have been obtained as = 3 and C = 4 0. Furthermore, 
we assume that the total bandwidth isW = 500 MHz, and the 
noise power spectrum density is Nq = —174 dBm/Hz. Denote 
by P the total transmission power, the average signal-to-noise 
ratio (SNR) at each receiving array element (without the lens 
applied yet) is then dehned as SNR= PE[/3]/(T^. We consider 
two communication environments characterized by different 
values of the maximum multi-path excessive delays T^, which 
correspond to: i) the narrow-band channel with Tm <C 1 /VF; 
and ii) the wide-hand channel with = 100 ns. 


power distribution 


with k}, = Ul 


In Fig. 0 the average spectrum efficiency is plotted against 
SNR for both the lens-based and the UPA-based mmWave 
MIMO systems in narrow-band communication, over 10^ 
random channel realizations. Note that for the UPA-based 
narrow-band MIMO system, the channel capacity is achieved 
by the well-known eigenmode transmission with WF power 
allocation based on singular value decomposition (SVD) over 
the MIMO channel matrix 1^. It is observed from Fig. 0 
that under the ideal AoA/AoD environment, the lens MIMO 
using OPDM achieves almost the same capacity as that by 
the conventional UPA-based MIMO. However, their required 
signal processing complexity and hardware cost are rather 
different, as will be shown in the next subsection. 


Fig. 0 compares the lens MIMO using OPDM versus 
the UPA-based MIMO using MIMO-OFDM in wide-hand 
communication. For MIMO-OFDM, the total bandwidth is 
divided into N = 512 orthogonal sub-bands, and a cyclic 
prehx (CP) of length 100 ns is assumed. It is observed in Fig.0 
that for the wide-hand communication case, the lens MIMO 



Fig. 6 : Capacity comparison of the lens MIMO using OPDM 
versus the UPA-based MIMO using eigenmode transmission 
in narrow-band mmWave communication. 



Fig. 7: Capacity comparison of the lens MIMO using OPDM 
versus the UPA-based MIMO using MIMO-OFDM in wide¬ 
band mmWave communication. 


achieves higher capacity than the UPA-based MIMO, which 
is mainly due to the time overhead saved for CP transmission. 

C. Complexity and Cost Comparison 

In this subsection, we compare the lens MIMO against the 
conventional UPA-based MIMO in mmWave communications 
in terms of signal processing complexity and hardware cost. 
The results are summarized in Table U and discussed in the 
following aspects: 

• MIMO processing: For the lens MIMO based mmWave 
communication, the capacity for both the narrow-band 
and wide-hand channels is achieved by the simple OPDM 
scheme, which can be efficiently implemented with signal 
processing complexity of 0{L), with O(-) representing 
the standard “big O” notation. In contrast, for the UPA- 
based mmWave MIMO communication, the capacity is 
achieved by the eigenmode transmission for narrow-band 














TABLE I; Complexity and cost comparison for lens MIMO versus UPA-based MIMO. 



Signal Processing Complexity 

Hardware Cost 

MIMO Processing 

Channel Estimation 

Antenna 

RF chain 

Narrow-band 

Wide-band 

Narrow-band 

Wide-band 

Lens MIMO 

om 

om 

OIL) 

0(^1 

M + Q 

2L 

UPA-based MIMO 

0{Mi}QvL) 

0(MvQvLN -b (Qu + Mv)N log N) 

O(MuQu) 

0{MvQvN) 

Mu + Qu 

Mu + Qu 


channel and approached by MIMO-OFDM for wide¬ 
band channel. The signal processing complexity for both 
schemes mainly arise from determining the eigen-space 
of the MIMO channel matrices, which has the complexity 
0 (MuQu niin{Mu, Qu}) for ^ generic matrix of size 
X Qu P^ - For a low-rank Mu x Qu channel 
matrix of rank L, the complexity can be reduced to 
0{M\jQ\jL) by exploiting its low-rank property |j^. 
Thus, the MIMO precoding/detection complexity for the 
UPA-based MIMO communication is 0{M\jQ]jL) and 
0{M\jQ\jLN) in narrow-band and wide-band communi¬ 
cations, respectively, where N denotes the total number 
of sub-carriers in MIMO-OFDM, which in general re¬ 
quires additional complexity of 0{{Q\j + M\j)N \ogN) 
at the transmitter and receiver for OFDM modula¬ 
tion/demodulation. As L ^ min{Mu,(5u} in mm Wave 
communications, the lens MIMO has a significantly lower 
signal processing complexity than the UPA-based MIMO, 
especially for the wide-band communication case. 

• Channel estimation: It follows from CD that the lens 
MIMO using OPDM only requires estimating L par¬ 
allel SISO channels for both narrow-band and wide¬ 
band communications, which has a complexity 0{L). In 
contrast, the conventional UPA-based MIMO in general 
requires estimating the MIMO channel of size Mu x Qu 
for narrow-band communication, and N different MIMO 
channels each of size Mu x Qu for wide-band commu¬ 
nication using MIMO-OFDM0 

• Hardware cost: The hardware cost for mmWave MIMO 
communications mainly depends on the required number 
of transmitting/receiving RF chains, which are composed 
of mixers, amplifiers, D/A or A/D converters, etc. For 
the lens MIMO system, it follows from CD that only L 
receiving/transmitting antennas located on the focusing 
points of the L multi-paths need to be selected to operate 
at one time; whereas all the remaining antennas can 
be deactivated. This thus helps to significantly reduce 
the number of RF chains required as compared to the 
conventional UPA-based MIMO, as shown in Table ||] in 
detail. 

IV. Lens MIMO under Arbitrary AoAs/AoDs 

In this section, we study the mmWave lens MIMO in 
the general channel with arbitrary AoAs/AoDs, i.e., the spa¬ 
tial frequencies {4‘R,u4'T,i}f=i are not necessarily integer 

^Note that by exploiting the channel sparsity in mmWave communications 
with small L, the channel estimation in UPA-based MIMO can be imple¬ 
mented with lower complexity via jointly estimating the multi-path param¬ 
eters \ai , (pud ; 4>t.i ; Tjl/Ui , which, however, requires more sophisticated 
techniques as in |21|. 


multiples of the spatial frequency resolutions of the receiv¬ 
ing/transmitting lens arrays. In this case, the power for each 
multi-path signal in general spreads across the entire antenna 
array with decaying power levels from the antenna closest to 
the corresponding focusing point. Let A > 0 be a positive 
integer with which it can be practically approximated that 
|sinc(a:)p « 0, V|x| > aQ It then follows from 0 and 0 
that the receive/transmit array responses for the Zth path are 
negligible at those antennas with a distance greater than A 
from the focusing point (see Fig. [^, i.e., 

aR,m{4>R,i) = v^^sinc(TO - « 0, Vm ^ Mi, 

(16) 

o-T,q{(l>T,i) = \/^sinc ((7 — Dt^t,i) ~ 0, Vg ^ Qi, 

(17) 

where Mi and Qi are referred to as the supporting receiv¬ 
ing/transmitting antenna subsets for the Zth path, which are 
defined as 

Mi = M ■.\m-bR^R^i\< , (18) 

Qz = {9 GQ:|7 -^t^t./|<a}, l = (19) 

Consequently, the (m, g)-th element of the channel impulse 
response matrix H(/) in (|^ has practically non-negligible 
power if and only if there exists at least one signal path I 
such that m G Ml and q G Qi- Since L ^ min{M, Q} due 
to the multi-path sparsity in mmWave systems, it follows that 
H(/) is in practice a (nearly) sparse matrix with block sparsity 
structure, where each non-zero block corresponds to one of 
the L multi-paths and has approximately 2A x 2A entries 
around the element {mi, qi), as illustrated in Fig. Note that 
depending on the AoA/AoD values, {Mi}fGi (or {Qi}fLi) 
may have non-empty intersection for different paths, i.e., 
certain antenna elements may receive/transmit non-negligible 
power from/to more than one signal paths, as illustrated by 
Q2 and Q 3 in Fig. 

Let Ms = Ml and Qs = iJiLi Qi be the supporting 
receiving/transmitting antenna subsets associated with all the 
L paths, and Hs(t) G be the sub-matrix of 

the channel impulse response H(/) corresponding to the 
receiving antennas in Ml 5 and transmitting antennas in Qs- By 
deactivating those antennas with negligible channel powers, 
the input-output relationship in ([^ then reduces to 

rMsib = Hs(t)®XQ,^(t) +ZMs(t) ( 20 ) 

L 

— X! ^l^R,Ms{4>R,l)^T,Qs{4>T,l)xQg{t - Ti) + Z,Ms 
1=1 

( 21 ) 

"^For practical applications, A = 1 is a reasonable choice, since 

|sinc(3i)p < 0.047, V|a)| > 1. 
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Fig. 8: An illustration of the sparsity for a lens MIMO 
channel with L = 3 paths. Dt = Du = 10. G 

{0.36,-0.27,0.08}, G {-0.2,0.12,0.24}. A = 1. 
Ml = {3,4},M 2 = {-3,-2},M 3 = {0,1},Qi = 
{-2}, Q 2 = {1, 2}, and Q 3 = {2, 3}. Note that Q 2 n Q 3 ^ 0 
due to the small AoD separation between path 2 and path 3. 


where G respectively denote the 

sub-vectors of r, and z in (|^ corresponding to the receiving 
antennas in Ms; and a'r Qg,XQg S denote the sub¬ 

vectors of ax and x corresponding to the transmitting antennas 
in Qs, respectively. 

Remark 1: It follows from that for mmWave lens 
MIMO system with arbitrary AoAs/AoDs, only |Ms| ^ M 
receiving and | Qs \ <C Q transmitting RF chains are generally 
needed to achieve the near-optimal performance of the full- 
MIMO system with all M -f Q antennas/RF chains in use. 
Furthermore, since |Ms| < SrLi \^i\ ~ 2AL, and |Qsl A 
Yld=i 12d ~ 2AL, the total number of RF chains required 
only depends on the number of multi-paths L, instead of the 
actually deployed antennas M and Q. 

A. Transceiver Design Based on RDM 

In this subsection, by exploiting the reduced-size channel 
matrix in we propose a low-complexity transceiver 

design based on PDM (instead of OPDM due to arbitrary 
AoAs/AoDs), which is applicable for both narrow-band and 
wide-hand mmWave communications. With PDM, L indepen¬ 
dent data streams are transmitted in general, each through 
one of the L multi-paths by transmit beamforming/precoding. 
Specifically, the discrete-time equivalent of the transmitted 
signal XQg (t) can be expressed as 

L 

^Qs W = 

1^1 

where n denotes the symbol index, si[n] ^ CAf (0,1) 
represents the i.i.d. CSCG distributed information-bearing 
symbols for data stream I, with transmit power pf, and 


^^T,asi4>T,i)si[n\, ( 22 ) 


ar, Qs ) / V At denotes the unit-norm per-path MRT 

beamforming vector towards the AoD (px^i of path 1. Note that 
we have used the identity jjaT’ .Qsi^T,lW ~ \\^t{4>T,i)V — 
At, V/. At the receiver side, the low-complexity per-stream 
based detection is used, where a receiving beamforming vector 
V; G with ||v/|| = 1 is applied over the receiving 

antennas in Mg for detecting s;[n]. Thus, we have 

sz[n] = vfr^g[n], ; = 1, • • • ,L, (23) 

where r 2 vig[n] is the discrete-time equivalent of the received 
signal r^g(f) shown in (|2T|. 

Next, we analyze the performance of the above proposed 
PDM scheme for wide-hand communications. The analysis 
for the special case of narrow-band communications can 
be obtained similarly and is thus omitted for brevity. For 
simplicity, we assume that the multi-path delays can be ap¬ 
proximated as integer multiples of the symbol interval Tg, i.e., 
Ti = niTg for some integer ni, VI. For notational conciseness, 
let ax,I = aT,Qs((/>T.i) and a^^i = an,Ms{^R,i)’ Based 
on and ( |2^ , the discrete-time equivalent received signal 
^MsIr] expressed as 

L 

^Ms W = X! '^k^R,k^T,k^Qs [r - Rk] + ZMs [r] (24) 



(25) 


Note that in we have decomposed the received signal 
^Ms[r] from the perspective of the /th data stream, which 
includes the desired signal component propagated via the /th 
path with symbol delay n;, the ISI from the same data stream 
received via all other L — 1 paths with different delays, and the 
inter-stream interference from the other L — 1 data streams. 
By applying the receiver beamforming in ( [2^ and treating 
the ISI and the inter-stream interference both as noise, the 
effective SNR for the /th data stream can be expressed as ( |2^ 
shown at the top of the next page. The achievable sum-rate 
is then given by i? = X]i^ilog 2 (l + 7;)- In the following, 
the two commonly used receiver beamforming schemes, i.e., 
MRC and MMSE beamforming, are studied to gain insights 
on the proposed PDM transmission scheme. 

1) MRC Receive Beamforming: With MRC, the receiver 
beamforming vector v; for data stream / is set to maximize 
the desired signal power from the /th path, i.e., = 

aR.i/s/Axi, V/. By substituting into ( |26| ), the SNR can 
be expressed as 0 shown at the top of the next page. Note 
that we have used the identity « Ar, V/. For two 

different paths /' /, define the transmitter and receiver side 
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7i = 


piAT\ai\^\vf^SLji^i 




yi. 


^mrc ^ 


Pi\ai\'^ ArAt 


Y^k^l ArAt 1'^'=!^ 


^R,yR,k 


^T,yT,l 


E Li 

I'^l Z^fe=i 


Pi' 

AjiAt 




^R,yR,k 


^T k^T,l' + O' 


(26) 

(27) 


inter-path contamination (IPC) coefficients as 


o“' - 

Pt — 


^T,yT,l' 


Al 


<1, Pr = 


^R,yR,I' 




< 1 . (28) 


The SNR in ( |27l l can then be simplified 


as 


MRC 

7i = _z. 


Pi\ai\ 


Efe^iPllQfclVflPT +EiVlEfe=lPl'|afePPHPT' + 

(29) 

__ Pi|Qi|^ ___ 

Efc^iPilafcPPflPT +E^^iPi' (lai'PPfl + I^iIVt ) + 

(30) 


where the approximation in ( |30l l is obtained by keeping only 
the two dominating inter-stream interference terms in ( |29] l with 
either k = I' or k = 1. 

It is observed from @ that for wide-band mmWave lens 
MIMO systems using PDM and the simple MRC receiver 
beamforming, the ISI is double attenuated as can be seen from 
the IPC coefficients and at both the transmitter and the 
receiver sides, and the inter-stream interference is attenuated 
through either transmitter-side IPC coefficient p“ or receiver- 
side IPC coefficient . Based on ([^, we have 


2 



aT,qi^T,l)aT,q{4'T,l') 

q^Qs 


2 


sinc{q - DT^T,i)sinc{q - Dt^t,I') , 

q&Qs 


(31) 

(32) 


which vanishes to zero for sufficiently separated AoDs such 
that 107’^ — (f>T,i'\ > 2A/_Dj’, or equivalently Qi n Qp = 0. 
Similarly this holds for the receiver side IPC coefficient . 
In Fig. M the IPC coefficient p“ is plotted against the AoD 
difference \(j)T,i — 4 >t,i' \ for different AoD resolutions provided 
by the transmitter lens array, which verifies that the IPC 
vanishes asymptotically with large AoD separations and/or 
high AoD resolutions. 

In the favorable propagation environment with both suffi¬ 
ciently separated AoAs and AoDs such that p^ « 0 and 
p“ « 0, yi' ^ I, both the ISI and the inter-stream interference 
in ( |29l l vanish. As a result, the SNR for the fih data stream 
reduces to 7 ; = pi\ai\‘^A rAt/ y, VZ, which is identical to 
that achieved by the OPDM in the ideal AoAs/AoDs case 
shown in Fig. In this case, PDM with simple MRC receive 
beamforming achieves the channel capacity for both narrow- 
band and wide-band mmWave communications. 



Fig. 9: Transmitter-side inter-path contamination coefficient 
versus the AoD difference \4>t,i — with different Dr in 
lens MIMO systems. 


2 ) MMSE Receive Beamforming: In the general case where 
the transmitter- and/or receiver-side IPC coefficients are non¬ 
zero due to the limited AoA/AoD separations and/or insuf¬ 
ficient AoA/AoD resolutions provided by the lens arrays, 
the PDM scheme suffers from both ISI and inter-stream 
interference, which needs to be further mitigated. One simple 
interference mitigation scheme is via MMSE beamforming at 
the receiver, for which the beamforming vector v/ in ( [ 2 ^ for 
the Zth data stream is set as EtI 


MMSE _ ^^R,l 


l=h 


^R,l 


(33) 


where C; is the covariance matrix of the effective noise vector. 
Based on ( [25] ), C; can be obtained as ( |35| ) shown at the top 
of the next page, where aR^k — aR^k/VA r, k = 1, • • • ,L. 
The corresponding SNR can be obtained as 

^MMSE yail'^piAra^aR,i (36) 

= \o^i\^Pi^R.i ('^Pi\o^k\^PTaR,kaR^k 

\ k^l 

+ aR^i 

l'^lk=l T R J 

(37) 

> 7 i“^S Z = l,--- ,L, (38) 

where ( [^ follows from the fact that x^A'^x > l/(x^Ax), 
V|jx|| = 1. The inequality in ( |38] l shows that MMSE beam- 






































II 


L L L 



L L L 

PiiA tA ji\ak\^ a_R,/ca-|j,fe + (T^I,V/, (35) 

I'^l k—1 


forming in general achieves better performance than MRC, 
since it strikes a balance between maximizing the desired 
signal power and minimizing the interference. In the favorable 
scenario with both sufficiently separated AoAs and AoDs such 
that « 0 and « 0 , VZ' 7 ^ I, it can be shown that the 
MMSE and MRC receive beamforming vectors are identical. 


B. Path Grouping 

As can be seen from ( [29l l and the performance of 
the RDM scheme with MRC or MMSE receive beamforming 
depends on the ISI and inter-stream interference power via the 
IPC coefficients and p*^, VZ' 7 ^ Z. In this subsection, the 
PDM scheme is further improved by applying the technique 
of path grouping, by which the paths that are signihcantly 
interfered with each other are grouped and jointly processed. 
It is shown that the PDM with path-grouping achieves the 
channel capacity for both narrow-band and wide-hand lens 
MIMO systems, provided that either the AoAs or AoDs (not 
necessarily both) are sufficiently separated. 

1) Sufficiently Separated AoAs: We first consider the case 
with sufficiently separated AoAs for all paths such that 
\^R,i — VZ' f Z, but with possibly close AoDs 

for certain paths, rfis may correspond to the uplink commu¬ 
nications where the receiving lens antenna array equipped at 
the base station has a large azimuth dimension (Dji ^ 1 ) 
and hence provides accurate AoA resolution; whereas the 
transmitting lens array at the mobile terminal can only provide 
moderate AoD resolution. In this case, it follows from ( fTS] ) that 
AdinAii' = 0, VZ' Z, i.e., form a disjoint partition 

for the supporting receiving antenna subset AAs- As a result, 
© can be decomposed into 

= CKia.ii^Mii4‘R,l)^T,Qsi4>T,l)^asit “ p) + 

l = (39) 


where G + ^ are respectively the sub¬ 
vectors of ^Ms in ( |2T] ) corresponding to 

the receiving antennas in Aii. ( [39l l shows that each receiving 
antenna only receives the signals via one of the multi¬ 
paths, thanks to the sufficient AoA separations such that 
the signals from different multi-paths are focused at non¬ 
overlapping receiving antenna subsets. However, the signal 
transmitted by certain transmitting antennas may propagate 
via more than one paths due to the possible overlapping of 
the supporting transmitting antenna subsets for different paths. 


Such a phenomenon is illustrated in Pig. 10 


Due to the single path received by each receiving antenna, 
the path delay t; can be compensated by the antennas in A4i. 



Pig. 10; An illustration of the effective channel in mmWave 
lens MIMO system with sufficiently separated AoAs. The path 
delays are labeled for each link. Gray antennas represent those 
with negligible power and hence can be deactivated. Note that 
path 2 and path 3 are grouped since they have similar AoDs 
at Tx (but different AoAs at Rx). 


As a result, ([39ll is equivalent to 


rMi' 

_ 

ai^RMi {4’R,l)^T,Qsi^T,l) 

XQs + 

ZAI 1 

^Ml. 


piL^R,ML i4’R.L)^T,Qs 


f‘ML_ 



In other words, with sufficiently separated AoAs, the original 
multi-path channel in ( |2T] i is essentially equivalent to a simple 
MIMO AWGN channel given in ( |40| , regardless of narrow- 
band or wide-hand communicationsljThe channel capacity of 
( |40| ) is known to be achieved by the eigenmode transmission 
with WE power allocation based on the MIMO channel matrix 
Hg. However, a closer look at H 5 reveals that it is still a 
sparse matrix due to the sparsity of the transmitting response 
vectors SLT^Qg((j)T,i), which can be further exploited to 
reduce the complexity for achieving the capacity of the MIMO 
channel in (|40li. 


^Recall from Section 


AoDs, i.e., 


IV-A 


: 0 or py 


that with either sufficiently separated AoAs or 
0, VP I, the ISI can be completely eliminated 


by PDM. 
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Fig. 11; An illustration of the effective channel in mm Wave 
lens MIMO system with sufficiently separated AoDs. The path 
delays are labeled for each link. Gray antennas represent those 
with negligible power and hence can be deactivated. Note that 
path 2 and path 3 are grouped since they have similar AoAs 
at Rx (but different AoDs at Tx). 


Recall that the transmitting array response vector 
aT,Qs{4>T,i) has essentially non-zero entries only for 
those transmitting antennas in the subset Qi C Qg. The 
main idea for the proposed design is called AoD-based path 
grouping, by which the L paths are partitioned into G < L 
groups such that paths I and I' belong to the same group if 
the transmitter-side IPC coefficient > 0, or equivalently 
if Qi n Qu ^ 0. Denote by Cg C {I,-- - ,L} the path 
indices in the 5 th group, p = 1, • • • , G. For instance, for the 
system shown in Fig. 10 we have G = 2 and Ci = {1} and 
£2 = {2,3}. In addition, denote by Qg = U/g£g Qi and 
Mg = 1Jzg£ Ml, g = 1, - • ■ ,G, the supporting transmitting 
and receiving antenna subsets for all paths in group g, 
respectively. By construction, {Qg}‘^^i and form 

disjoint partitions for the supporting transmitting antenna 
subsets Qs and Mg, respectively. Therefore, the input-output 
relationship in ( |40l i can be decomposed into G parallel MIMO 
AWGN channels as 


^Mg — -I- , g — I," ■ , G, (41) 

where and Xg^ S denote 

the sub-vectors of r^g, z^vig and xgg in ( |40| ), respectively; 
and Hg = Y.i^Cg^i^R,Mg{(t>R,MT,Qg^^T,i) denotes the 
corresponding MIMO channel matrix for group g. The channel 
capacity of •ED is then achieved by the eigenmode trans¬ 
mission over each of the G parallel MIMO channels, which 
have smaller dimension and hence require lower complexity 
as compared to the original channel in ( |40l l without path 
grouping. 

2) Sufficiently Separated AoDs: Next, we consider the 
case with sufficiently separated AoDs for all paths such that 


|0T,z — ^T,i'\ > 2/SIDt, or Qi n Qv — f I', but with 

possibly close AoAs for certain paths. This may correspond 
to the downlink transmission with accurate AoD resolution 
(Dt ^ 1 ) at the base station transmitter, but with only 
moderate AoA resolution at the mobile terminal receiver. In 
this case, {Qi}ffi form a disjoint partition for the transmitting 
antenna subset Qg, and the input-output relationship in (| 2 T| 
can be re-written as 

L 

= X! ^l^R,Ms{4>R,MT,Qi{4’T,lMQi{t - T;) -f 
1=1 

(42) 


The expression in ( |42| l shows that the signals sent by each 
transmitting antenna arrive at the receiver only via one of the 
multi-paths, as illustrated in Fig. 11 This thus provides the 
opportunity for path delay pre-compensation at the transmitter 
by setting the transmitted signal as xg, (f) = X g^ {t + r/), VI. 
As a result, ( |42l i can be equivalently written as EU shown at 
the top of the next page. 

Similar to the previous subsection, ED shows that with 
sufficiently separated AoDs, the lens MIMO system is equiv¬ 
alent to a I Ads I X |Qs| MIMO AWGN channel. This holds 
regardless of narrow-band or wide-hand communications. The 
channel capacity of ( |4^ is achievable by eigenmode trans¬ 
mission with WF power allocation based on the equivalent 
channel matrix Hs. Similar to Section IV-B1 1 by exploiting 
the channel sparsity of Hs, we can design a low-complexity 
capacity-achieving scheme by employing the AoA-based path¬ 
grouping at the receiver side. Specifically, the L signal paths 
are partitioned into G < L groups such that paths I and I' 
belong to the same group if their supporting receiving antenna 
subsets have non-empty overlapping, i.e., A4i n Me f 0 - 
Denote by £g C {1, •••,£}, g = 1,---,G, the subset 
containing all paths in group g. For instance, for the system 
shown in Fig. 11 we have G = 2 and £1 = {1} and 
£2 = { 2 ,3}. In addition, denote by Qg = IJzg£ 2z and 
Mg = Uz 6 £ Ml, 5 = 1, • • • , G, the supporting transmitting 
and receiving antenna subsets for all paths in group g, respec¬ 
tively. Similar to Section IV-Bl the input-output relationship 
in ( |4^ can then be decomposed into G parallel MIMO AWGN 
channels as 


^Mg = Hgx'g^ + ZAi,, P = 1, • • • , G, (44) 

where G clAtglxi g (j^|Cg|xi (jgjjote 

the sub-vectors of r^g, z^g and ^Cs “ ED’ respectively; 
and Hg = Y.i^Cg‘^i^RMgi.^R,MT,Qg^^T,i) denotes the 
corresponding MIMO channel matrix for group g. The channel 
capacity of ( |44| is then achieved by the eigenmode transmis¬ 
sion over each of the G parallel MIMO channels each with 
reduced size. 


C. Numerical Results 

In this subsection, we evaluate the performance of the 
proposed PDM in a wide-hand mmWave lens MIMO system 
by numerical examples. We assume that the lens apertures at 
the transmitter and receiver are At = 100 and An = 50, 
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rA^s = ■■■ Cy.L^R,Ms{4’R,L)^T,QS^T,L)\ 


Hs 


+Z^ 


^QiJ 

‘ss 


(43) 


respectively, and the azimuth lens dimensions are Dt — 
20 and = 10, respectively. Accordingly, the number 
of transmitting and receiving antennas in the lens MIMO 
systems are Q = 41 and M = 21, respectively. For the 
benchmark MIMO system based on the conventional UPAs, 
the number of transmitting and receiving antennas are set as 
Q\] = 400 and My = 200, respectively, for achieving the same 
antenna apertures as the lens MIMO system. For both the lens 
MIMO and UPA-based MIMO systems, antenna selections are 
applied by assuming that the number of RF chains at the 
transmitter and receiver are Mrf = Qrf = 2AL, where L 
is the number of multi-paths and A is a design parameter 
to achieve a reasonable balance between performance and 
RF chain cost. We set A = 1 in this example. For the 
lens MIMO system, the AoA/AoD based antenna selection 
method given in ( [T8] l and ( [T9| ) are applied at the receiver and 
transmitter, respectively. However, since the optimal antenna 
scheme for the UPA-based MIMO-OFDM system is unknown 
in general, we adopt the power-based antenna selection due 
to its simplicity and good performance | [25| . We assume that 
the mmWave channel has L = 3 paths with AoDs given 
by (j)T,i G {—15°, 10°,45°}, which are sufficiently separated 
based on the criterion specified in Section 1V-B2| On the other 
hand, the AoAs of the L paths are assumed to be equally 
spaced in the interval [—$fj/2, ^ii/2], with referred to as 
the AoA angular spread. Furthermore, the maximum multi- 
path delay is assumed to be = 100 ns and the total 
available bandwidth is IF = 500 MHz, which is divided into 
N = 512 sub-carriers for the UPA-based MIMO-OFDM. The 
CP length for the OFDM is set as 100 ns. 

In Fig. [T^ the spectrum efficiency achieved by different 
schemes is shown for the mm Wave communication with AoA 
angular spread = 150°. Note that for simplicity the 
power allocation for the PDM with MRC and MMSE 

receive beamforming is obtained via WF by assuming L 
parallel SISO channels with power gains {|ai| It 

is observed from Fig. 12 that the UPA-based MIMO-OFDM 
gives rather poor performance, which is expected due to the 
limited array gain with the small number of antennas selected. 
In contrast, the lens MIMO systems with the three proposed 
PDM schemes achieve significant rate improvement over the 
UPA-based MIMO-OFDM with the same number of RF chains 
used or antennas selected. Moreover, Fig. 12 shows that in 
the low-SNR regime, PDM with the simple MMSE and MRC 
receive beamforming achieves the same performance as that 
with path grouping, which is expected due to the negligible 
inter-path interference in the low-SNR regime. While as the 
SNR increases, the three PDM schemes show more different 
performances due to their different interference mitigation 
capabilities. The performance gaps are more pronounced for 



Pig. 12: Average spectrum efficiency achieved by various 
schemes in wide-hand mmWave MIMO communication with 
antenna selection. The number of transmitting/receiving RF 
chains are Mrf = Qrf = 6. AoA angular spread is = 
150°. 



Fig. 13: Average spectrum efficiency achieved by various 
schemes in wide-hand mmWave MIMO communication with 
antenna selection. The number of transmitting/receiving RF 
chains are Mrf = Qrf = 6. AoA angular spread is $/£ = 10°. 


systems with smaller AoA separations, as shown in Fig. 13 
for = 10° as compared to Fig. 12 for = 150°. 
This implies the necessity of more sophisticated interference 
mitigation techniques (such as path grouping) for PDM when 
the paths are severely coupled with each other due to the 
limited AoA/AoD separations. 
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V. Conclusion and Future Work 

In this paper, we studied the use of lens antenna arrays 
for mmWave MIMO communications. The array response 
of the lens antenna array was derived and compared with 
that of conventional UFA without the lens. We showed that 
the proposed lens antenna array significantly reduces the 
signal processing complexity and RF chain cost as compared 
to the conventional UFA in mmWave MIMO communica¬ 
tions, and yet without notable performance degradation. We 
proposed a new low-complexity MIMO spatial multiplexing 
technique called FDM, for both narrow-band and wide-hand 
communications. Analytical results showed that the FDM 
scheme is able to achieve perfect ISI rejection as long as 
the AoAs or AoDs are sufficiently separated, thanks to the 
energy focusing capability of the lens antenna. Finally, for 
cases with insufficient AoA/AoD separations, a simple path 
grouping technique was proposed for FDM to mitigate inter¬ 
path interference more effectively. 

There are a number of interesting directions that are worthy 
of future investigation, which are briefly discussed as follows. 

• Elevation AoAs/AoDs: For systems with non-negligible 
elevation AoAs/AoDs, the array configuration of the lens 
antenna arrays needs to be refined. In this case, the 
antenna elements should be generally placed on the focal 
surface of the EM lens to exploit the elevation angular 
dimension as well, instead of on the focal arc only as 
considered in this paper. As a result, the signal multi¬ 
paths can be further differentiated with the additional 
elevation AoA/AoD dimension. 

• Multi-User Systems: The FDM for the point-to-point 
mmWave MIMO communication can be extended to the 
general path division multiple access (FDMA) for multi¬ 
user mmWave systems, by which a number of users with 
well separated AoAs/AoDs can be simultaneously served 
with low-complexity and low-cost transceiver designs. 
The transmission scheduling of users based on their 
AoAs/AoDs is also worth investigating. 

• Channel Estimation: In this paper, we assume per¬ 
fect channel state information at both the transmitter 
and receiver, while in practical mmWave systems such 
knowledge needs to be efficiently obtained via well- 
designed channel training/estimation/feedback schemes. 
For mmWave MIMO communications with conventional 
arrays, channel estimation is a challenging task due to the 
large-antenna dimension as well as the low SNR before 
beamforming is applied Q l[T§, pT[ ; whereas with 
lens antenna arrays, by exploiting its energy focusing 
as well as the multi-path sparsity of mmWave channels, 
the effective channel dimension is significantly reduced 
and the pre-beamforming SNR is greatly enhanced. 
Therefore, channel knowledge can be obtained far more 
efficiently as compared to conventional arrays, which 
deserves further study. 


Appendix A 
Froof of Lemma[T] 

To derive the array response of the proposed lens antenna 
array given in Lemma [T] we first present the fundamental prin¬ 



Fig. 14; Top view of a planar EM lens placed in the y-z plane 
with focal point Bo{F, 0, 0) for normal incident plane waves. 


ciple of operation for EM lenses. EM lenses are fundamentally 
similar to optical lenses, which are able to alter the propagation 
directions of the EM rays to achieve energy focusing or beam 
collimation. 

Eig. 14 shows a planar EM lens of size Dy x Dz placed 
in the y-z plane and centered at the origin. Denote by Bq 
with coordinate {F, 0,0) the focal point of the lens for normal 
incident plane waves, where F is known as the focal length. 
The main mechanism to achieve energy focusing at Bq is to 
design the phase shift profile which represents the 

phase delay provided by the spatial phase shifters (SFS) of 
the lens at any point (0, y, z) on the lens’s aperture, such that 
all rays with normal incidence arrive at Bq with identical phase 
for constructive superposition |35|. We thus have 


^{y,z)+kod{y,z,Bo) = fl>o, 


yiy,z) e 


Dy 

Dy 

X 

Dz 

Dz' 

2 

2 

2 

2 


(45) 


where ko = 2 tt/X is the free-space wave number of the 
incident wave, with A denoting the free-space wavelength, 
d{y,z,Bo) = y/W -f 2 /^ -f is the distance between the 
point (0, y, z) on the lens’s aperture and the focal point Bq, 
and $0 is a positive constant denoting the common phase delay 
from the lens’s input aperture to the focal point Bq. The phase 
shift profile is then designed to be 


$(?/, z) = $0 - kQs/F'^ -f -hz2, 


V(y,z) G 


Dy 

Dy 

X 

Dz 

Dz 

2 

2 

2 

2 


(46) 


As can be seen from ( |46] l, due to the different propagation 
distances from the lens’s aperture to Bq, the phase shift profile 
varies across the lens apertures with different y and z values. 
In general, larger phase delay needs to be provided by the SFS 
located in the center of the lens than those on the edge. 

With the phase shift profile designed in ( |46l l to achieve focal 
point Bq for normal incident wave, the resulting phase delay 
from the lens’s input aperture {0,y,z) to an arbitrary point 
B{xB,yB,ZB) is then given by 


i^iy,z,B) = <^{y,z)kQd{y,z,B), (47) 


























15 



It is interesting to observe from ( |5^ that with the spatial 
phase shifting provided by the EM lens, the resulting signal at 
the focal arc of the lens can be represented as the Fourier 
transform of the arriving signal s{y) at the lens’s input 
aperture, with 0 G [—1,1] and y G 


-Dl2,Dj2 


(|^ 


referred to as the spatial frequency and the spatial time, 
respectively. 

For uniform incident plane waves with azimuth AoA 
or equivalently with spatial frequency f = sin(0), as shown 


in Fig. 


15 


we have s{y) = 


equivalenfry 




-xo{(j))e 


.j^ysinicj)) 


Ky) = 


y/ffyU, 




J27r2/ 


(54) 


where xo{(j>) is the input signal arriving at the lens center 
with AoA (j), and y/DyD^ is a normalization factor to ensure 
Fig. 15; A planar EM lens placed in the y-z plane with oblique that the total power captured by the lens is proportional to its 
incident plane wave of azimuth AoA f. effective aperture A = DyDz/X^- By substituting ( |54l l into 

(|53l), we have 


where d{y, z, B) = \Jx\G- {ys - vY + {zb - z)"^ denotes 
the distance from the point (0, y, z) on the lens to point B. 
Of particular interest is the held distribution on the focal 
arc of the lens, which is dehned as the arc on the x-y 
plane with a distance F from the lens center, as shown in 
Fig. 15 Let B{Fcos9,—Fsm6,0) be a point on the focal 
arc parameterized by angle 9 G With ( |46l l and 

we have 


r{9) = xo{(jf)'/Asinc (^D{9 — , 0e[—1,1]. (55) 

It then follows from ( |55] l that the effective lens response on 
its focal arc for incident plane waves with AoA f (or spatial 
frequency (p) is 

ag{(j)) ='/Asmc{D{0 — f)), 9 g[—1,V\. (56) 


f{y,z,9) 


4-0 - + y2 + z^ + ko^/Fi 


4-0 + koysmO, 


+ j /2 -(- ^2 _|_ 2yF sin 9 

(48) 

(49) 


For the lens antenna array with the mth element located at 
position Bm{Fcos{9m), —Fsm{9m,),Q), it follows from ( |56| ) 
that the array response can be expressed as 


where ( |49| ) follows from the hrst-order Taylor approximation 
and the assumption that F ^ Dy,Dz. 

Let s{y, z) denote the incident signal arriving at the lens’s 
input aperture. Due to the linear superposition principle, the 
resulting signal on the focal arc of the lens can then be 
expressed as 



where in ( |5T] i, we have assumed that s{y, z) = s{y), V(j/, z) G 
_i^ D^ X [—11^, 11^], which is true for uniform incident 
plane waves with negligible elevation AoAs. For notational 
convenience, we assume that $o = 2n7r for some integer n, 
so that it can be ignored in ( [ST] ). Furthermore, by dehning 

D=^, y=j, 9 = sin{9), (52) 

the relationship in ( [5T] l can be equivalently written as 

m^Dzl ^ ~s{y)e-^'^^~^ydy, 0G[-1,1], (53) 

J-D/2 


where s{y) with y G 


-DI2,DI2 


is a linear scaling of the 


arriving signal s{y) given by s(y) = As(Ay). 


= VAsinc \^D{sin9m — sin(/))j , Vm. (57) 

In particular, with the critical antenna spacing specihed in 
0’ the array response in ( |57| ) reduces to 

amif) = "x/Asinc (jn — Dsinf^ , Vm. (58) 
This completes the proof of Lemma 
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